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Abstract 

In this paper we describe a systematic approach for constructing asymptotic boundary 
conditions for isotropic wave-like equations using local operators. The conditions take a recur- 
sive form with increasing order of accuracy. In three dimensions the recursion terminates and 
the resulting conditions are exact for solutions which are described by finite combinations of 
angular spherical harmonics. First, we develop the expansion for the two-dimensional wave 
equation and construct a sequence of easily implementable boundary conditions. We show 
that in three dimensions the analogous conditions are again easily implementable in addition 
to being exact. Also, we provide extensions of these ideas to hyperbolic systems. Namely, 
Maxwell’s equations for TM waves are used to demonstrate the construction. Finally, we 
provide n um erical examples to demonstrate the effectiveness of these conditions for a model 
problem governed by the wave equation. 
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1 Introduction 


The topic of boundary conditions for wave-like equations has been an active area for the 
past two decades, with much work focusing on the development of sequences of conditions 
of increasing accuracy based on local operators. Growing interest in demanding problems in 
computational electromagnetics and aeroacoustics has increased the need for more accurate and 
efficient techniques. In this paper we discuss and derive a sequence of boundary conditions 
for isotropic problems. Preliminary results of this work are found in [8]. Generalizations 
to anisotropic problems, such as the linearized compressible Euler system, will be discussed 
elsewhere. The present development originates with the well-known reference [2]. It differs 
from that and from others in that conditions of arbitrary order are expressed recursively using 
auxiliary functions. This allows them to be implemented with great ease. It also allows the 
adaptive determination of the order of the boundary condition, which we will demonstrate 
elsewhere [7]. Moreover, in three space dimensions, the recursion truncates for finite spherical 
harmonic expansions, leading to an exact condition. This condition is similar to that of Grote 
and Keller [3, 4], but is somewhat easier and cheaper to use as it avoids spherical harmonic 
transformations. 

This paper is divided into four parts. In section 2 we describe the procedure for the 
two-dimensional wave equation. The key idea is to construct a sequence of operators that 
approximately annihilate the residual of the preceding element in the sequence, viewed as a 
function on the artificial boundary. The sequence begins with the first order operator proposed 
in [2]. The recursive operators are constructed in such a way that they yield a stable class 
of conditions. The stability issue however is not discussed in this paper and will be reported 
elsewhere. Also, we indicate how these conditions can be approximated numerically. In section 
3 we discuss analogous constructions for the three-dimensional wave equation. Generalizations 
to a hyperbolic system, that is Maxwell’s equations in the TM case, axe presented in section 
4. Finally, in section 5, we give some numerical results. 

2 The Two-Dimensional Wave Equation 

A formal series representation for an outgoing solution of the two-dimensional 

1 d 2 u _ d 2 u 1 du 1 d 2 u 

<? dt 2 dr 2 r dr + r 2 dO 2 

takes the form: 

oo 

u = ]Ca n ( 0 )/ n (r,t), 

n=0 
oo 

/n = £ r _ fc _ i / 2/ „( rt _ r)) 
k - 0 


wave equation 
( 2 . 0 . 1 ) 

( 2 . 0 . 2 ) 

(2.0.3) 
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where, 

tin \/ + 2^ 2 n2 fn n _ _ n 2 Q 

(/fc+i) 2(/t + 1) ■'*” n 

(2.0.4) 

By setting 

1 du du 1 

1 1 U = VJ \ , 

eft 9r 2r 

(2.0.5) 

it is easily observed that 

OO 



U>1 = °nW/ n,1 ( r > t ). 

n— 0 

(2.0.6) 


oo 

Z”’ 1 = £ r ~ k ~ 3/2 'Tk,ifk(ct - r), 
k = 1 

(2.0.7) 


7fc,l = - fc - 

(2.0.8) 


By direct computation we note that uq = 0(r~ 2 )u = 0(r 5/2 ). We now seek relations 
approximately satisfied by w\, which in addition do not involve radial derivatives. In [8] we 
proposed a sequence of time derivatives of w\ to derive higher order conditions. While they are 
accurate, we found them to be only marginally stable. Here we propose adding a damping term 
to the time derivative of wi, as first suggested by Barry, Bielak and MacCamy [1]. Precisely, 
we apply the operator ^ to w\. This leads to the following calculations: 


c at r ^ 


n=0 k=0 

oo oo 


OO OO 1 1 

= E E + b 2 - " 2 ) - 

n=0 k= 0 

1 1 d^u i V — fc-5/2 (\(i> _1_ __ ^ k __ /n 

= 2^302 +£ a "z] r 8^’ 


n=0 k— 1 

The remainder term is further reduced if S\ is chosen to eliminate the first term: 

CO 

$ 1 = 1 , V>2=Y1 a nf n ’ 2 i 

n = 0 


^ j 

f n ' 2 = r- k - 5/2 i k , 2 fk(ct - r), 7fc,2 = ~ *)• 


(2.0.9) 

( 2 . 0 . 10 ) 


k = 2 


Clearly, we formally have w 2 = 0(r 4 )u = 0(r 9//2 ). 



We now state and prove the general formula. Let 


w i = Y^o, n (d)f n ’ 3 (r,t), 


71—0 


f nJ = k j 1/2 7k,jfk(ct - r), 

k=j 


where 


Then, for j > 2: 


7ktJ = (-iy2 l -ik(k-l)---(k-(j-l)). 


1 dwj j 

i 4- = 


u - i) 




+ ~ w i = - I „2 f — ^-1 + 


1 , . _ n 

c dt ' r~ J 4r 2 ' 4r 2 dB* +“ , J+ 1 » J “ A- • • 

To prove this we simply expand every term. On the left-hand side we have: 

00 fi. 4. 1\2 _ „2 oo 

- £ r-»4-VW, + £ r -‘-Wj-KjfS- 


h=j~l 

On the right-hand side we have: 


k=j 


oo (A _ IV _ r) 2 oo 


k—j — l " k=j+l 

Separately we check k = j — 1, k = j and k > j 4- 1. 


= (-!>• 


2 j 


F-Ui-J 

2j 


= (— l) J-1 2 2 ~-*(j — 1)!— — 2 i 2 - ■ -- 


— Tj-lj-l' 


(J - ^) 2 - " 2 


0 + 5) 2 -” 2 , • 

2j + 2 ' +J7j ' j 


= (-l) J - I 2 1 - J (j + 1)! ° + 2 j - ^ 2 — ~ 


= 7W-1- 


U - h) 2 - n 2 


( 2 . 0 . 11 ) 

( 2 . 0 . 12 ) 

(2.0.13) 

(2.0.14) 

(2.0.15) 

(2.0.16) 
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3 The Three-Dimensional Wave Equation 


The wave equation in three dimensions is given by: 


1 d 2 u d 2 u ^ 2 du ^ J_ V7 2„ 

c 2 dt 2 dr 2 r dr r 2 s ’ 

(3.0.21) 

where the spherical Laplacian is given by: 


1 d ( . n du\ , 1 d 2 u 

v,u ~ smeae \ sm) a») sm 2 ed<t?' 

(3.0.22) 

We seek a formal solution of the form: 


CO n 

«=£Ean,m(M)/ n ’ m M), 

n— 0 m— 0 

(3.0.23) 

On,m = (An ,m COS m<£ + Bn,™ sinm<£)P™(cOS 0), 

(3.0.24) 
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(3.0.25) 


r m ( r,l) = £r- t - 1 /""(cf-r). 
k = 0 

Assuming each term in the first summation to be a solution, and recalling the fact that 

we derive the following recursion relation: 

/ rn y m\t (k(k -hi) — Tl(ri j~ 1 )) . n ,m 

Uk+i) - 2(k 4- 1) 5 

where is arbitrary. Note that we may choose f£' m = 0 for k > n so that: 

n 

/ n - m (r, <) = £ r - fc -i/ fc n ’ m (ct - r). (3.0.28) 

fc= o 


(3.0.26) 

(3.0.27) 


Moreover, given f n ’ m (R,t), the right-hand side of (3.0.28) can be rewritten as an nth order 
differential operator on /” ,m , 


/ n ' m (i?,o = X> 

A;=0 


-k-l 




jn—k /n,m 


(3.0.29) 


r)n = 1, 


%- n 

j=fc+i 


2j 

n(n + 1) - j(j - 1) ’ 


k < n. 


(3.0.30) 


Hence, given f n ’ m (R,t), the expansion (3.0.28) certainly exists. We emphasize that we do 
not need to compute the functions f£’ m in order to implement the boundary conditions. (See 
(3.0.47)-(3.0.48).) 

We now repeat the derivation of the preceding section. The main difference will be the 
termination of the sequence for each finite n. Set, 


1 du du 1 

(3.0.31) 

7T7 + *3“ + ~ U = WU 
c at or r 

oo n 


Wl = £ £ «n,m(M)/ n ’ m> 1 (M), 

(3.0.32) 

n— 1 m= 0 


/“■"■•I = f; r -‘- 2 7M/?' m (c <- r) . 

(3.0.33) 


fc= l 


By direct computation we note that w\ = 0(r 2 )u = 0(r 3 ) and 


7*1,1 = -*• 


(3.0.34) 
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Now consider: 


1 dw\ 
c dt 


w 1 


r 


LXJ ft 't 

-EE “».™ E i--*' s ((* + !)(/?«)' + 

n=l m=0 fc=0 

oo n n -i 

E E <Vm E r-^HUkik + 1 ) - »(n + D) - ( 3 . 0 . 35 ) 


n=l m=0 fc= 0 

1 do n 


V 2 ,u + E E E + 1 - Ml )/» 


rn 5 m 


n=l m=0 


fc=l 


The remainder term is further reduced if 6\ is chosen to eliminate the first term: 


$1 = 1, W 2 = £ £ °n,m/' 


n,m,2 


n=2 m= 0 


(3.0.36) 


yn,m,2 = £ r -*-3 ~ r), 7*, 2 = ^(* “ 1)- (3-0-37) 

fc =2 

Clearly, 102 = 0(r _4 )u = 0(r -5 ). Moreover, if we were in the special case of a solution in the 
span of the spherical harmonics of index less than 2, then w 2 = 0 and the condition is exact. 
We now consider the general case. Suppose for 1 < j < p, p > 2, 

00 n 

Wj = £ £ (3-0-38) 

n=j 7n=0 


f n ' m ' 3 = £ r- k ~ 3 - 1 ') k,jfk’ m (ct - r), 

k=j 


where again 7 fc 5 j is as given in (2.0.13). Then: 


and 


1 df n 


,77l,p 


c <9t 


_). P fnjnjp _ y r 


k—p 

= £ 

Jt=p-1 


7fc,p^(/jt ) +-/* J 

(■ 


Tfe+jjg 
2(k -t- 1) 


(3.0.39) 


= — (;nT + ;' w ) • < m40) 

° ^ n=pm=0 N 


(3.0.41) 


(k(k + 1) — «(n + 1)) + P7fc,p ) fk ’ , 
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(3.0.42) 


where we have used the fact that 7 p -i, p = 0. Noting that: 


7fc+i, P 
2 (k + 1) 


1 

— ^7fe )P - li 


and 


we find 


1 dw p 
c dt 


, P 

+ -w p 
r 


(k-(p- 1)) 

7 k,p — 0 7&,p-li 


(3.0.43) 


^2 ( V * + P(P~ (3.0.44) 

oc n n 

+ E E^™ E r- fc -P- 2 ^-(fc-p)(fc-(p-l))/ fc n ' m . 


As the k = p and k = p — 1 terms in the final summation vanish, the n — p — 1 and n — p 
terms are eliminated. Moreover, 


^~^(k -p){k - [p - 1)) = 7t lP +i, 


so the remainder is given by: 


oo n 


w P + 1 = E E °"- m E r k v 2 7fc )P +i fk ,m - 

n— p-f 1 m— 0 A:=p-fl 


By induction we obtain: 


1 du du 1 

+ - + -u = wi, 
c ot or r 


1 dwj j 1 

+ J -Wi = — 


c dt 


4r 2 


(v 2 4- j{j - 1)) Wj-i + Wj+i, j = 2,... 


where u>o = 2 u, Wj satisfies (3.0.38)-(3.0.39) and, hence, 

Wj = 0(r -2j )it — 0(r -2j,_1 ). 


(3.0.45) 

(3.0.46) 

(3.0.47) 

(3.0.48) 

(3.0.49) 


Using p auxiliary functions involves setting to p +i = 0 which leads to an expected accuracy of 
0(r~ 2p ~ 2 ). However, in the special case: 

f n ' m = 0, n > p, (3.0.50) 

that is when the solution is a finite sum of spherical harmonics, we have: 


w p+ i = 0 . 


(3.0.51) 
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Hence, in this case, the boundary condition defined by p auxiliary functions is exact. Therefore 
it is an exact, local boundary condition in space and time in the same sense that the conditions 
proposed by Grote and Keller [3, 4] are exact and temporally local. The necessity of a spherical 
harmonic expansion in the Grote-Keller formulation probably renders theirs less efficient than 
the formulation proposed here. We note that a distinct formulation with the possibility of 
avoiding the spherical harmonic expansion is mentioned in their work, but is not implemented. 


4 The Two-Dimensional Maxwell System 


We now consider the generalization of this technique to a hyperbolic system. Namely, we study 
the two-dimensional Maxwell system in the TM case, given by: 


d_ 

dt 




E z 

H x 

\Hy 


+ >>;! 



(4.0.52) 


where 

/ 0 — £ sin 6 j cos 6 \ 

R= — i sin 9 0 0,0 

V i cos0 0 0 / 


0 — jcos# — jsin# ^ 

— i cos 6 0 0 

— ^sin# 0 0 j 


(4.0.53) 


Note that it is possible to directly apply results for the wave equation to this system as in 
[5]. However, it is instructive to proceed directly. Setting c = (e/i) -1 / 2 , the progressive wave 
expansion now takes the form: 


f E z \ 00 

H x = £ r~( k+1 Mf k (ct - r )a k (0), (4.0.54) 

\Hy J *=0 


where 

f' k = h-u (4-0-55) 

{cl + R)a 0 = 0, (4.0.56) 

{cl + R)a k — —{k — l/2)Ra k -i + (4.0.57) 

It is easily verified that the matrix on the left-hand side of (4.0.56) is singular, so that a 
nontrivial expansion exists. Generally one must solve an eikonal equation to determine the 
form of the expansion. A convenient orthonormal basis for expanding the solution is provided 

° ^ 

— sin# . (4.0.58) 

cos# J 


oy: 


ei = 


1 \ 

0 

w 


p = 



T = 


9 



These satisfy the relations: 


Re i = — r, Rp = 0, Rt 

i 

= -ei, 

(4.0.59) 


e 


_ 1^1 

0ei = p, ©p = — e 1? 

p e 

©T = 0, 

(4.0.60) 

II 

II 

i 


(4.0.61) 


In terms of this basis, a right nullvector, u, and a left nullvector, i, of (cl + R) are given by: 


u — e i 




The solution of (4.0.56) is: 


ao = A(9)u , 

where A(0) is an arbitrary function. For (4.0.57) we expand 


and note: 

The equations become: 


a* — Eke i + 



Eq = — Co = A, Bo — 0. 


E k + C k = - {(k - 3/2)C fc _x + B£_ x ) , 


(4.0.62) 

(4.0.63) 

(4.0.64) 

(4.0.65) 

(4.0.66) 


E k + C k = -(k- 1/2 )£*_!, (4.0.67) 

B* = (4.0.68) 

Clearly, the system is solvable if and only if 

{(k - 3/2 )C*_i + = (k - l/2)E k ~i, (4.0.69) 

and then it only determines E k + C k , B k . The complete determination of E k , C k follows 
from the imposition of (4.0.69) at the next order. (Note that (4.0.69) does hold when k = 1.) 
Solving these leads to the simple recursion relations for A: > 1: 

E k = -(2 k)~ l (E'U + (k- 1 /2) 2 £*_x) , 

B/: = —E k _h 
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(4.0.70) 

(4.0.71) 



Ck = (2k)- 1 (E'U - ( k 2 - 1/4)jE?*_i) • 
To determine a boundary condition we first note that: 

(4.0.72) 

OO 

E z = 'Z r ~ {M/2) fk( ct ~ r)Ek, 

fc= 0 

(4.0.73) 

I OO 

H r = cos 9H X + sin 9H y = J- E r~( k+1 ^ f k (ct - r)#*, 

’ ^ fc=0 

(4.0.74) 

j OO 

H$ = — sinflifx + cos 9H y — J- E r~^ +1 ^/jfc(ct — r)Cfc. 

’ ^ fc=o 

(4.0.75) 

To leading order we use (4.0.65) to obtain: 


j OO 

= E ^" {fc+1/2) /fc(c< - r)(JS? fc + C fc ) = 0(r -1 ) • |(£ z , fl*, H y )\. 

(4.0.76) 


fc=l 


A more accurate formula follows from taking the time derivative and using the recursions: 


E z +J-H e ) = Y, r ~ {k+l/2) fk(ct-r){E k +i+C k +i) = — E r~( k+1/2) f k (ct-r)(k+l/2)E k , 

CCft V 6 fc =0 r k = 0 

(4.0.77) 

which implies 

~(E, + ^H e ) + ^E z = 0(r- 2 ) ■ \(E z ,H x ,H y )\. (4.0.78) 

To find higher order conditions, we must derive a relationship for the evolution of the 
right-hand side, w \ , of the equation above. For j > 1 set 


Wj = kjE k = 0(r~ 2 i) ■ \(E z ,H x ,Hy)\. (4.0.79) 

k=j 


Here ikj is again given by (2.0.13). Note that uq is indeed the right-hand side of (4.0.78). 
Differentiating Wj with respect to t and using (4.0. 55), (4.0.70) yields: 


k=j 

00 

= E r ~ k ~ j ~ Z ^ 2 fklk+ijEk+i (4.0.80) 

fc=j- 1 

1 OO 

- r- k -i-V 2 f k 'r k+ i, j ( k +V- 1 (EZ + (k + l/2) 2 Ek). 

1 k=j—\ 
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From the identity 


(k + 1/2) 2 = (k- ( j - 1 ))(k - j) + 2 j(k - (j - 1)) + (j - 1/2) 2 , 


and (2.0.13) we finally derive 


1 dwj 
c dt 


1 

4 r 2 


( g 2 Wj- 1 

^ d9 2 


+ U - 1 / 2 ) 



— -Wj + Wj+l . 


(4.0.81) 


(4.0.82) 


Putting these together we have a boundary condition whose formal order of accuracy is 


r —2m. 


1 dWj 
C dt 


+ ~Wj = 


\jt {Ez + ^i He) + b Ez ~ wu 

1 ( d 2 wj-i 


4r 2 { d0 2 


+ (j - l/2) j! Wj- 1 J + Wj+ 1, j = 1, . . . ,m - 1, 

wo — 2 E z , = 0 . 


(4.0.83) 

(4.0.84) 

(4.0.85) 


Of course in numerical implementations, this condition must be supplemented with some 
extrapolation of the “outgoing” variable, E z - y^iftan and a computation of the variable 
characteristic along the boundary, ^normal- (For a circular boundary these are H tan = Hg , 
Ff normal = H r .) A three-dimensional version also exists. Moreover, it is exact in the same way 
that the condition described in the previous section is exact. See [7] for more details. 


5 Numerical Experiments 


We consider numerical experiments with these conditions in the case of the two-dimensional 
wave equation with c = 1. Further calculations along with experiments for Maxwell’s equations 
are in progress [7]. The test case considered here is a Dirichlet problem exterior to a disk of 
uni t, radius. A range of tangential Fourier modes is excited by a time-dependent amplitude 
function which both oscillates and decays. Precisely we take: 

u(r = 1, 9, t) = 1 cosn 0’ (5.0.86) 


for n = 0, 5, 10, 25, 50. The modal wave equation governing the amplitude of cos nO, that is: 


d 2 u d 2 u 1 du n 2 

dt 2 dr 2 ^ r dr r 2 U ’ 


(5.0.87) 


is discretized by a standard second order central difference scheme in space-time and the 
boundary conditions given by equations (2.0.17)-(2.0.18) are discretized as described below. 
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To construct a second order approximation to these equations, we apply the boundary 
condition half way between the last two mesh circles and we also use the fact that w J+ 1 is 
small compared to Wj, (w\ small compared to u). Then we approximate (2.0.17) by: 


1 

2 St 


(( U N + u N-l) ~ ( u iV + u N-l)) 


+ + U N) ~ \ U N-1 + 

+ + «V-1 + ut N + u N-\) 

3 t 1 t - 
= - 2™ 1 ’ 


- A / / 


and (2.0.18) by 


~ wt j) + + w ^ 


2 R 

0-1/2 ) 2 -n 2 


SR 2 


{Wj - 1 + 1 ) 


+ 


-w 


1 t- 


'i+i 


— -w 


'j+ 1- 


For wo we use the sum of u at N and N- 1. Note that although the system would formally be 
implicit in the ^-derivative approximation for a non-modal algorithm, no matrix factorizations 
axe involved if we solve sequentially. That is, x is computed before its ^-difference is used to 
update Wj. The calculations below me for t = 15 with varying artificial boundary location R. 
The value of m is the number of auxiliary functions used, so that the formal asymptotic order 
is 2m + 2. The value of 6r is chosen to yield about four digits of accuracy in the numerical 
solution. We are thus forced to use a much finer mesh for n = 50 than for n = 0, with the 
number of points increasing by a factor of 30. The ‘exact’ solution is computed with the 
same mesh spacing but R = 18. In all cases we choose 8t = 0.8<5r. Errors were measured at 
t — . 5 , 1 , 1.5, . . . , 15. The time station at which the maximum error occurred, not included in 
the tables, was always t = 15 for n = 0, but earlier in the other cases. The largest error was, 
however, always recorded after the waves with largest amplitude passed through the artificial 
boundary. We note that at later times the solution is fairly small, so that the relative L 2 error 
criterion used here is quite strict. Its precise definition is: 


2 /.\ _ X)j( u long( r j ; t) u i r ji 0) 2 

Zj(uion g (rj,t)) 2 


(5.0.88) 


We first consider n = 0 where 5r = 5 x 10 3 . 
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Maximum Relative Error n = 0 


m 

R = 2 

CO 

II 

R = 4 

0 

3.8 x 10" 1 

5.0 x 10' 1 

5.4 x 10" 1 

2 

2.2 x 10 _1 

1.2 x 10- 1 

4.0 x 10" 2 

5 

2.9 x 10 -i 

3.1 x 10~ 3 

1.7 x 10~ 4 

10 

4.4 x 10~ 4 

1.2 x 10 -4 

1.7 x 10~ 4 

20 

4.6 x 10~ 5 

1.2 x 10" 4 

1.7 x 10~ 4 


The second case is given by n = 5 with dr = 2.5 x 10 3 . 


Maximum Relative Error n = 5 


m 

R = 2 

i? = 3 

-Sf 

II 

05 

0 

5.0 x 10" 1 

4.5 x 10 _1 

4.0 x 10" 1 

2 

2.3 x 10“ 2 

6.0 x 10“ 3 

1.8 x 10" 3 

5 

3.9 x 10" 5 

6.6 x 10~ 5 

1.3 x 10" 4 

10 

3.9 x 10“ 5 

6.6 x 10~ 5 

1.3 x 10~ 4 

20 

3.9 x 10~ 5 

6.6 x 10- & 

1.3 x 10“ 4 


The third case is given by n = 10 with dr = 1.6 x 10 3 . 


Maximum Relative Error n = 10 



R = 2 

R = 3 

R = 4 

0 


1.1 

1.1 

2 

3.0 x 10~ 2 


1.5 x 10~ 3 

5 

5.7 x 10" 4 

mma 

QQom 


3.9 x 10" 5 

OSQQ&I 



3.9 x 10" 5 

9.6 X 10~ 5 

2.0 x 10" 4 


The fourth case is given by n = 25 with dr = 3.3 x 10 4 . 














Maximum Relative Error n = 25 


m 

R = 2 

R = 3 

R = 4 

0 

5.4 x 10 _i 

5.7 x 10 _1 

5.7 x 10" 1 

2 

4.5 x 10 -3 

7.0 x 10 -4 

1.9 x 10" 4 

5 

2.1 x 10~ 5 

6.1 x 10“* 

1.2 x 10~ 4 

10 

2.1 x 10" 5 

6.1 x KT 5 ' 

1.2 x 10~ 4 

20 

2.1 x 10 -5 

6.1 x 10" 5 

1.2 x 10~ 4 


Finally we take n = 50 with Sr = 1.6 x 10 4 . 


Maximum Relative Error n = 50 


m 

R = 2 

R = 3 

R = 4 

0 

4.1 x 10 _1 

4.5 x nr 1 

4.5 x 10 _1 

2 

2.3 x 10“ 3 

3.7 x 10~ 4 

1.3 x 10~ 4 

5 

2.2 x 10" s 

6.6 x 10 _;> 

1.3 x 10" 4 

10 

2.2 x 10 -5 

6.6 x 10" 5 

1.3 x 10~ 4 

20 

2.2 x 10' 5 

6.6 x 10" 5 

1.3 x 10“ 4 


With increasing m we rapidly drive the error down to the level of the discretization error, 
so that the results become m-independent. Note that this final error is not zero as we are 
comparing with a numerical long domain solution. Indeed, the final error (as m — > oo) may 
increase with increasing domain size. We further checked the m-independence of the results by 
carrying out additional experiments with m = 30, 40, 50, 60, which produced errors identical 
to the m = 20 case shown below. Note that for none of the domains tested were the m = 0 
results acceptable. For m = 20, in contrast, the errors are dominated by discretization errors 
even on the smallest domain. Generally, the n = 0 case is the most difficult, as predicted by 
the analysis of [9]. 

6 Concluding Remarks 

We believe the boundary conditions formulated in this paper provide an easily implementable 
path to high accuracy, at a computational cost which is favorable in comparison with currently 
available alternatives. There are, however, some open issues to be addressed. From a theoret- 
ical perspective, asymptotic error estimates are as yet unproven in the two-dimensional case, 
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and one as always expects some time nonuniformity in the approximation. From a practical 
perspective, our current formulation is restricted to circular /spherical boundaries. It would, of 
course, be useful to have more flexibility. The expansions themselves should be valid outside 
some circle/sphere containing sources, scatterers and other perturbations. However, the direct 
imposition of the conditions on a boundary of different shape requires normal derivatives of 
the auxiliary functions. It would then be necessary to define and store them off the boundary, 
which we want to avoid for reasons of cost and complexity. Generalizations to anisotropic 
problems such as the convective wave equation and the linearized compressible Euler equation 
are also possible. Now the boundary must take a special shape determined by the eikonal 
equation associated with the underlying velocity field. This will be discussed in detail in a 
later paper. 
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